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Abstract 

The present review is devoted to the problems of finite-size scaling due to the 
presence of long-range interaction decaying at large distance as l/r d+a , where d is 
the spatial dimension and the long-range parameter a > 0. Classical and quantum 
systems are considered. 



1. Finite-size systems and critical phenomena 

A common wisdom is that the singularities in the thermodynamic functions at a critical 
point may occur only in the thermodynamic limit. If the system is fully finite (as any real 
system is), no such singularities exist and, strictly speaking, no phase transitions occur. 
These facts naturally give rise to the following questions: 

a) Why the bulk theory turns out to be adequate to the experimental evidence for 
finite objects? 

b) How the singularities appear when no phase transitions occur in any fully finite 

system? 

The answer to the first question lies in the fact that a macroscopic body is close, under 
some circumstances, to the idealized thermodynamic limit. 

The answer to the second question becomes evident if one recalls that the limit function 
of a sequence of analytic functions needs not to be analytic. Therefore, what happens is 
that in the thermodynamic limit some thermodynamic functions of the system develop 
singularities which are attributes of phase transitions. So, in a finite system one expects 
an appreciable rounding of the critical point singularities. 

The bulk correlation length £oo{T) measures the distance over which particles in a 
system are significantly correlated. Then, as the temperature T approaches the critical 
temperature T c , £oo(^) diverges. When £,oo(T) attains a magnitude of the order of the 
characteristic size L of the finite system, then the boundary particles at the opposite sides 
of the system become well correlated, and ordering cannot continue to build up further 
in the restricted dimensions. It is certainly reasonable to expect that the rounding of the 
phase transition is controlled by the criterion: £oo(^) ~ L. 

The description of this rounding and crossover was first formulated by M. Fisher (1972) 
and is called Finite-Size Scaling (FSS)( see e.g. P). 

An important application of FSS is to analyze numerical data obtained from simu- 
lations on relatively small finite systems, and, thereby to obtain knowledge of the bulk 
systems. 

Finite-size effects are much stronger in systems of particles with long-range (LR) 
interactions because every particle directly feels the influence of the boundaries. 
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To study finite systems, one must at the start address two basic issues, namely the 
overall geometry of the system and the specific nature of the boundary conditions. 

Geometry: We will use the general notation L d ~ d x oo d and distinguish three impor- 
tant cases (d,d' are spatial dimensions): d! — - a fully finite system, d! — 1 - a system 
with the geometry of a cylinder, and d' — (d — 1) - a system with slab (or film) geometry. 

Boundary conditions: Here and bellow we will consider periodic boundary conditions 
(p.b.c). Close to the bulk critical point for a d-dimensional system with a finite size L 
one has for the free energy density 

/«(D = /(r) + o(e-^). (1) 

where f(T) = foo(T) is the bulk free energy density. In this case there are no surfaces, 
edges and corners. 

The experimental observability of FSS effects is hampered by the fact that nature does 
not provide us systems with p.b.c. The rapid exponential approach to the thermodynamic 
limit is one of the reasons why p.b.c. are preferred in Monte Carlo simulations. One 
observes that the l.h.s. in the above expansions is a regular function of T whereas the 
r.h.s. is not: the bulk density /(T), for example, is singular at the critical point T c . 
Therefore, Eq.(JIJ can hold only away from T c . Obviously, for T ~ T c one needs an 
alternative, (FSS) formulation. 

2. Finite-size scaling hypothesis 

Scaling hypotheses are made in the form of statements about properties of thermodynamic 
quantities in terms of homogeneous functions. 

To simplify the further discussion we consider that: the studied system is below its 
upper critical dimension d u , the infinite system has a second order phase transition at 
a critical temperature T c , and the system has length L in all finite directions. The 
violation of these restrictions will cause complications and modifications of the standard 
FSS hypothesis pQ. 

For a physical quantity P with singular behavior in the thermodynamic limit the 
content of the standard FSS hypothesis is to assume the existence of a scaling function 
P L (T) such that: 

Pl(T) ~ Poo(T)F £ (L/e(r)). (2) 

As far as L is finite, Pl(T) must be a regular function of T. The (universal) function 
F L (x) must compensate the singularity of P 00 (x). Eq.Q may be written in equivalent 
form: 

P L (t) ~ If-'TLitL 1 /"), (3) 

where p x - the critical exponent of the observable Pl, v - the critical exponent of the 
correlation length and t = (T — T c )/T c . The hypothesis is expected to be valid when L is 
large and T approaches T c . 

What is the common status of the above hypothesis? One can usually meet the 
following statements, e.g. J. Cardy, in 0: "FSS is by now well-established theoretically, 
at least for systems with short-range interactions or quite recently Chen and Dohm 
in [3| : "The finite-size scaling hypothesis asserts ... in the absence of LR interaction ..." 
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We would like to point out the restrictions related with the nature of the interaction - 
one insists that it must be short-range (SR). In this situation the problem is to find out: 
does the FSS take place if the interaction is of the LR type? As we will explain bellow 
there is a positive answer to this question. 

We will consider interactions which decay algebraically at large distances with the 
standard notation 

u(r) ~ -l/r d+<T , a > 0. (4) 

For a < the interactions are nonintegrable, and so, under standard definitions, the 
thermodynamic limit does not exist. The case of nonintegrable interaction is beyond the 
scope of the present study. 

The LR interaction (J3J) enters the expressions of the theory only through its Fourier 
transform pQ. The corresponding small - q expansion of the Fourier transform has the 
general form 

= v + v 2 c[ 2 + v a \q\ a + w(q) < a ^ 2, (5) 

with u>(q)/|q| CT — > 0, for q — > 0, i.e. in the long-wavelength approximation. Further on 
we will formally relate a = 2, to the SR interaction since then (jSJ) is the Fourier transform 
of an interaction decaying exponentially with distance. 

Depending on whether in (JHJ) a is less or bigger than 2 we will speak about leading or 
subleading LR term respectively. 

Nota bene: The subleading LR term does not affect the bulk critical behavior, but 
influences the finite-size one (for details see |H 03 E|)- 

The defined above LR interaction has to be a mimicry of the real one. Nevertheless, 
from the pure theoretical point of view the considered LR interaction seems to be of 
specific interest. The reason is that renormalization group (RG) predictions obtained 
on the basis of the e-expansion, can be verified on ideal testing ground. The upper and 
lower critical dimensions are 2a and a respectively. Since a is a continuous parameter 
the value of e = 2o~ — d or e = d — a would be small enough for integer values of the 
dimensionality. This places us in the rare situation when the outcome of the computer 
simulations, obtained by means of the Monte Carlo method, can be directly compared to 
the predictions of the e-expansion with e arbitrarily small. 

But along this line of consideration there are some problems. 

Firstly, although the general scheme of the e-expansion has been widely accepted, the 
part of it related to LR systems has not been well understood until now (see [7j and 
references therein). 

Secondly, the study of LR systems leads to prohibitively large computational require- 
ments. Still recently a novel Monte Carlo algorithm which has an efficiency that is inde- 
pendent of the number of interactions per spin has been announced jHJ |§] . 

Thirdly, Monte Carlo analyzes of the RG predictions typically apply FSS concepts. 
However, specific problems arise in the proper generalization of the FSS concepts in the 
case under consideration jTUl UJ] • 

3. Problems with the correct definitions of the correlation length 
By definition the correlation length is 

£i(T) = - lim[i2/lnG? 00 (R;T)]. (6) 
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Alternatively, one may consider the second moment of the bulk pair correlation function 
and define the effective correlation radius: 



These two most commonly used definitions of the bulk correlation length are unam- 
biguous and equivalent (up to a constant) in the case of exponential with the distance 
R — |R| decay of the bulk pair correlation function l7oo(R; T). A distinctive feature of the 
LR interaction is that the function C7oo(R; T) decays as R~ d+cr when R — > oo for T > T c , 
and both definitions yield £i(T) = £2 CO = 00 if a < 2. 

Anyhow complications with correlations that decay as power laws in a whole domain 
of thermodynamic parameters, rather than only at special points, arise even in the bulk 
case, provoking the idea of "generic scale invariance" in the theory of phase transitions 
(see, e.g. The FSS under consideration encounters the same phenomena mixed 

with specific boundary effects. One can overcome the difficulty with the absence of good 
definition of the correlation length (when a < 2) following the ideas proposed in ^UJ > 
where instead of £(T) a bulk characteristic length A(T) is used. This characteristic length 
determines the length-scale of variation of the correlation function and diverges at the 
critical point. 

A finite system has three characteristic length scales: a linear size L, a length A^(T) 
which determines the scale of variation of the correlations, and a microscopic length a 
(e.g., a is the lattice spacing). The approach proposed in [TT] involves the following steps: 

First step: It consists of three assumptions. 

The first one defines the finite-size characteristic length Xl{T) from the large- distance 
asymptotic behavior of the finite-size pair correlation function Gl(R, T). 

The second assumption is the homogeneity of any thermodynamic quantity of the 
finite-size system, such that its bulk limit is singular at T = T c , as a function of the two 
dimensionless ratios A^(T)/a and L/a. 

The third assumption is the existence of finite thermodynamic limits for the charac- 
teristic length and the pair correlation function. 

Second step: It concerns the relationship between the finite-size length Xl{T) and its 
bulk limit Aoo(T); the corresponding homogeneity assumption involves the dimensionless 
ratios Aqo (T)/a and L/a. 

4. Mathematical problems 

Classical case: In the case of SR interaction the following replacement is used as an 
indispensable part of the FSS calculations (see, e.g. |T]). 



where q and q are d-dimensional and one-dimensional discrete vectors, respectively. This 
is the so called Schwinger parametric representation. The aim of this replacement is two- 
fold: to reduce the <i-dimensional sum to the corresponding effective one- dimensional one, 
and to give the dimensionality d the status of a continuous variable. 




nl/2 
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In the case of leading q a term in (jSJ), one cannot use the Schwinger representation in 
its familiar form. Then the following generalization of (jSJ) has been suggested [13] : 



E 



i 



m z + \q\ 



m 



dtQ a (m 4/a t) 



1 d 



(9) 



where the function Q a (t) for < a < 2 is given by 



Qa{t) 







dyexp(-ty)Q a (y), Q a {y) 



sin(<j7r/2)|/ 



a/2 



^From (jUJ) and (fTU]) one obtains 

1 



E 



m 2 + |q| CT 



dtQ a 
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7T 1 + 2y CT / 2 COs(cT7r/2) + ?/° 
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tm 2 + |q| 



(10) 
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The identity (JTTJ) demonstrates the possibility to reduce the LR case to the SR case 
with an integration over an additional parameter. So, all the mathematical machinery 
developed for the SR case, in principle, may be used. Let us note however that such 
nonlocal expressions are unsuitable for explicit calculations. It is more convenient to use 
the relation |14j : 



Qa{x) 



a/2,a/2 



(-W 2 ) 



12) 



where E a ^(z) is an entire function of the Mittag-Leffler type defined by the power series 



Eaj>(z) = E 



' T(ak + (3) 



a > 0. 



(13) 



and to study the finite-size properties of the system resulting from the analytical properties 
of E a3 . 

Pure quantum case: At T = 0, the quantum counterpart of © is 



E 



i 



X [m 2 + |q 



all/2 



7T JO 



dpV — ^ 7, — r~ j— • 

q m z + p A + |q| CT 



(14) 



Equation (]T4*]) displays the well known property of the pure quantum case. The aux- 
iliary variable p 2 acts effectively as an extra dimension. Indeed the pure quantum system 
corresponds to a d+1 dimensional classical system with the geometry of a cylinder L^xoo 1 . 
Introducing m 2 (p) := m 2 + p 2 with the help of Eq. (jll|) from Eq. (]T4*]) one gets 



E 



q [m 2 + |q 
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(15) 



The integral representation ()15j) illustrates the mathematical difficulties appearing in the 
pure quantum case. It is shown that by two additional integrations the problem can be 
effectively reduced to the classical SR case. 

On the other hand, the following modification of (jHJ) has been prosed ^5] 



E 



q 



all/2 



I ' dtC a {m 4/(T t) 



i d 



Y,exp(-q 2 t) 



(16) 



where 

C a (x) = x^G^/^-x^ 2 ), (17) 
i.e. the function C a (x) is related to the entire function defined by the power series 

1 00 V(k 4- 1 IT\7 k 

G " w ^Sw a>0 -' 3>0 ' (18) 

Some analytical properties of G a /2,a/i{z) (see also [IB]) may be established with the 
help of the identity 

Ga/2,a/4:{z) = ~ J Q E ^/2,a/2 {-{z + p 2 )) dp, (19) 

obtained from Eq.(JT3J) and the relation of its l.h.s and r.h.s. with the functions Gv/2,a74 (2) 
and E (T /2 t a-/2(z), respectively. 



5. Theoretical models 

Here we will mention the models and the related works in which leading LR interactions 
are considered in the context of FSS. 

The mean spherical model was considered in ^UJEIEIIE], |17j-|22j. The Hamiltonian 
of the model is 

H = ~J2j ij S i S j -HY,Si, (20) 

ij i 

where Si is the spin variable at site i, Jij is the interaction matrix between spins at 
sites i and j, and H is an ordering external magnetic field. The spins obey the spherical 
constraint J2i(<~>i) — N, where (• • •) denotes the standard thermodynamic average taken 
with the Hamiltonian TC and N is the the total number of spins. 

The quantum version of model (|2(J|) was considered in [22] • The Hamiltonian of the 
model is 

« = | - \ E^i ^ + ^E^ 2 -#E<^ (2i) 

% ij i i 

where Si are spin operators at site i, the operators V% are "conjugated" momenta and 
[iSj,iS-] = iSa', [Pi,Vi'} = i5u' (with % = 1). Model in the large-n limit is equivalent 
to the quantum 0(n) nonlinear sigma model. 

The ci-dimensional 0(n)-symmetric model was considered in |21]-j2Hl- The model is 
defined by 

,2 „ 1 



d x 



V 3 ] + ?W + 2 M o^ 



22 



where if is a shorthand notation for the space dependent n-component field <p(x), ro = 
roc + to (t oc T — T c ) and uq are model constants. The first term denotes k a \ip{k)\ 2 
in the momentum representation. A strongly anisotropic version of (|2*2*|) with exponents 
depending on the direction was considered in j29j . 
The quantum version of model (|2*2*|) is 



U M = \ /q 1/T dT Iv ^ [( Va/ ^) 2 + + l UoV 



,2 „ 1 



(23) 



where p is a shorthand notation for the space-time dependent n-component field (p(x,r) 
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6. Verification of FSS 



Classical phase transitions 

Analytically the case a < d < 2a has been investigated first exactly in the spherical limit 

(n = oo) HH [13 HH, dl-llll and by the RG methods for n = 1 |25|, n > 1 |2E1E|. 

Finite-size effects are predicted in both regions of first order [T3 HH1 EH E] and second 
order phase transitions [TSJ UHl UU UH ESI EE1 Elj • As a result in the latter case it has 
been shown that the FSS relations can be written in form equivalent to Eq. (|3J) . Some 
thermodynamic quantities, for instance, the shift of the critical temperature due to the 
finite-size effects [2T], the susceptibility [2011111, the pair correlation function [TUHH] and 
the Binder's cumulant at the bulk critical temperature T c j2Hl EHl Elj and above it [2E1 E2j 
as a function of \fe up to C(e 3 ' 2 ) have been obtained. 

Away from the critical region, the critical behavior of the systems is dominated by its 
bulk critical behavior 

Nota bene: A distinctive feature of the LR case is that for tL l / v ^> 1 the finite-size 
corrections are not exponentially small as in the SR case; they vary instead as algebraic 
power of the variable tL x l v [El CHI EE! • 

Some problems exist with the comparison of RG analytical results and numerical data. 
It has been shown (see also [SO]) using Monte Carlo simulation (for n = 1), that the 
Binder's cumulant ratio is linear in e. The analytical evaluation of the Binder's cumulant 
for the 0(n) symmetric if 4 model, however, showed that it is linear in y/e. A possible way 
to resolve this controversy between the Monte Carlo method and the analytical results 
is to carry on finite-size calculations to higher loop order [23] • This could ameliorate the 
analytical results, which would be comparable to those obtained by numerical simulations. 
Here we would like to mention that higher loop corrections that are dealt with through 
minimal subtraction scheme (see jSUESEEl) and e-expansion have not been performed, 
even for the more simple case of SR interaction and therefore much knowledge is still 
absent. 

Quantum phase transitions 

Quantum phase transitions [33] occur at zero temperature as a function of some non- 
thermal control parameter such as composition or pressure and are driven by quantum 
fluctuations. In addition to being of an interest for various low-temperature phenom- 
ena, quantum phase transitions are important because they are believed to play a crucial 
role in quantum information science, e.g. in the phenomenon known as entanglement - 
the resource that enables quantum computation and communications (see [35] and refs. 
therein) . 

The main feature of quantum systems is the coupling between statics and dynamics. 
As a result the dynamical critical exponent z plays an important role in scaling of the 
"imaginary time dimension". For nonzero temperature this extra "dimension" extends 
only over a finite interval L T . Here L T ~ 1/T is "finite-size" in the temporal (imaginary 
time) direction and so the temperature's influence in the quantum critical point my be 
determined by FSS effects [1.. 

Unlike classical models, where scaling can be done uniformly for all "spatial" di- 
mensions, quantum models are anisotropic in general, and therefore the "space" and 
"imaginary-time" directions will not scale in the same fashion. According to the general 
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hypothesis of the FSS theory extended here to quantum (anisotropic) systems, a physical 
quantity A (r, L, T), which is singular in the thermodynamic (bulk) limit at the quantum 
critical point (r = 0), will scale like 



A(r, T, L) = FA S (rb 1/v , Tb z , bL' 1 ) . 



(24) 



Here, r measures the distance from the quantum critical point. In the scaling form (}2"4*j) 
p corresponds to the engineering dimension d + z of the system in the case when the 
scaling function refers to the singular part of the free energy. For the other physical 
quantities of interest it is the critical exponent measuring the divergence of the bulk 
thermodynamic function A at the critical point divided by v (for the correlation length 
p = 1, for the susceptibility p = j/v, etc). Depending on the choice of the RG rescaling 
factor b we obtain different scaling functions A s , which are related among each other by 
some appropriate change of the scaling variables [To] . 

The fact that the inverse temperature can be used as an additional size in the imaginary- 
time direction creates strong anisotropy in the system and in general, one needs a shape- 
dependent formulation of the FSS ^01I2H]- This will lead to the establishment of some 
changes in the scaling properties of the finite quantum system. In this case we can con- 
sider the quantum to classical and the finite-size to the bulk system [TH], and as well as 
other crossover phenomena. 

Results on the Casimir effect 

The correlation function serves as a measure of how fluctuations at one point are correlated 
with fluctuations at another point. The confinement conditions, imposed on a system with 
correlation function decaying as power law in space induce a LR force between the surfaces 
limiting the system. One can generally call this phenomenon Casimir effect. There can be 
different mechanisms leading to this phenomenon that are related with the scale invariance 
or mentioned above "generic scale invariance" . Prominent examples are systems at critical 
points or systems with spontaneously broken global continuous symmetry that leads to 
massless modes: " spine waves" or Goldstone bosons (statistical-mechanical Casimir effect 
[32J), and the discussed for the first time by H.G.B. Casimir (1948) constrained zero-point 
vacuum fluctuations of the electromagnetic field (quantum - mechanical Casimir effect ). 

In a system with geometry L 1 x oo^" 1 at or bellow the critical point a LR "statistical- 
mechanical Casimir force" [TJ EI21 EH] appears. According to the standard FSS hypothesis, 
near the critical point of the bulk system one expects for this force (the magnetic field 



where X(x) is an universal function and k is a nonuniversal metric factor. In the case of 
SR interaction the universal scaling structure of Eq. (J23J) has been confirmed by RG and 
exact calculations [TJ • 

The investigation of systems with LR interaction possesses some peculiarities: due 
to the character of the interaction there exists a natural attraction between the surfaces 
bounding the system. In this case only results within the spherical model with p.b.c. are 
obtained [22]- It is possible to formulate some general statements: the Casimir force is 
always negative, (i.e. it is a force of attraction between the surfaces bounding the system) 
for any T and a > 1, as well as for any T > T c if a < 1; the behavior of the Casimir 



H=0) 




(25) 
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force depends strongly on the range of interaction a, being (at T = T c ) a monotonically 
increasing function of a; in the neighborhood of T = T c , the Casimir force is an increasing 
function of T if a > 1, and possesses a complex behavior if a < 1. 

For the Casimir force there is an exact expression obtained in the framework of the 
quantum spherical model (see pQ) 

i 71 Casimir(T, A; L) ~ X Casimir (x , p), (26) 

with scaling variable x = L 1 ' v (A -1 — A^ 1 ) , and p = L Z /L T ,L T = A/T. Here A is the 
quantum parameter that governs the transition near the quantum critical point A c . 

The obtained exact expression for the universal scaling function Acasimir^, p) gives the 
possibility of analysis including issues as the sign of the Casimir force and its monotonicity 
as a function of the temperature and parameter a [231. 

The author thanks J. Brankov, H. Chamati and D. Danchev for valuable discussions. 
Support by Bulgarian NSF, Grant F-1402, is acknowledged. 
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